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A treatment based on the analysis of data from Satellite 1963 38c 
has been developed to calculate the decay lifetime f of the artificial 
electron belt produced by the Starfish nuclear explosion. The treatment 
makes use of experimental data and describes the decay life time t in 
the form of a three-dimensional continuum in B, L and energy space. The 
results are valid over the entire inner zone of the trapped radiation belt. 
Reference maps of constant t contours in B-L, B-E, and L-E space have 
been constructed and are discussed* 


Introduction 


Seven artificial radiation belts have been made by the explosion of 
hich- altitude nuclear devices since 1958. Of theoe the Starfish event of 
July 9, 1962, a low Intitule nuclear explosion, produced an intense region 
of trapped artificial electrons with very long lifetimes. 

The spatial and temporal characteristics of the artificial radiation 
zone have received considerable treatment in the literature, and the decay 
lifetimes for the injected electrons have been discussed by ocvcral authors 
r von Allen . 1964; Cladls et ol ., 1965; Bostrom and Williams , 1965; Brown . 

:9 66; Mcllvaln . 1966; Ibullkas et ol .. 1967; Beall et al .. 1967; Pfltzer , 
1968; Bostrom et al ., 1970]. Bie picture that emerces is that of a fact 
decay at very low L values (L < 1*2), which is consistent with the prediction 
of atmospheric scattering theory f Walt , 1964]. Short lifetimes also occur at 
the outer edge of the inner radiation zone in the slot region. For the 
intermediate L values, the decay is slower with longest lifetimes occurring 
at L ~ 1.4. At low L values (1.15 < L < 1.4), a B dependence is also 
observed r Beall ct al .. 1967]* 

During the years since 1962, computer programs were developed at 
Goddard Space Flight Center to estimate the flux decay of artificially 
injected electrons and the radiation hazards to which satellites would be 
exposed r Stossinopoulos , 1965# 1967]* The decay lifetimes used in the 
calculation of orbit 'integrated, vehicle- encountered particle fluxes were 
based on fragmentary published data, and as such they .were necessarily 
restricted and not very accurate F Stossinopoulos , 1965]* In this paper 
we present a unified treatment for the decay of the artificial electrons 
of the inner zone by expressing the decay lifetime t as a continuous 
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function of I), L and energy E. The data uned in the analyolo In ccncntially 
taken from a circle satellite, 1963 38c, and covcro the time span from launch, 

September 28, 1963# to December 1968 vhen, It 5? felt, the energy channels of 
the catcllito measured ecoentlally the radiation of the natural backcround 
rBo ntrorn ct ol ., 1970]. The catellitc description and instrumentation can 
be found in the literature f Williams and Smith, 1965 ; Bostrom ct al ., 1967; 

Beall et al ., 1967 ]. 

Discussion 

In order to facilitate the handlinc of an otherwise difficult and 
complex problem, the following simplify ing assumptions ore made: 

(a) That t remains constant in time. 

(b) That solar cycle, diurnal, and secular- change effects can 
be neglected. 

(c) That the decay of the Starfish electrons can be approximated 
by an exponential form. 

(d) That the B-dcpcndencc of t docs not extend beyond L 1.1*0. 
Assumption (c) is Justified by the 1963 3^0 data; assumption Id) is supported 
by comparison with 0G0 data which measures a wider range of equatorial pitch 
angles than 38c [Pfitzcr, 1968; Pfttzer end Wincklor . 1968]. 

When ordering the data in terms of the variable L, three distinct regions 
emerge (see Fig. l), each identified by a characteristic functional dependence 
of t on the parameters B, L, and E. We denote the three regions us I, 

II, and III. In region I, which extends from L = 1.19 bo L « 1.1*0 (for 

all electron energies), *r is a function of all the variables: 

the field strength li, the magnetic shell parameter L, and the energy E. 

The B dependence of the decay constant has been previously discussed by 
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Beall ct nl ., (1967), and it appears to become Indeterminate for L ;» l.H. 
Consequently, in our treatment, ve will take t in region to depend only 
on L and E. Region II extends from L > l.fe to approximately L~ 2.3, with 
the latter boundary being energy dependent. Curves, obtained for region II 
from the APL data by plotting t vs E for all three channels at given 
Ir values, suggest that the high energy electrons exhibit the most rapid 
decay and that at the time of the 1963 3&C measurements almost no artifi- 
cials with energies E > 2 Mev remained in regions of L > 1.8. This is in 
general agreement with currently prevailing opinion P Pfltzer . 1968; Vettc, 
private communication, 1970]. Thus, in generating t vb E plots for 
region II (2.3>L>l.U)ve have not made use of the data from; the 1963 
38c channel 3 (E 2.4 Mev) and have Instead extended the dependence 
linearly to the higher energies, using only data obtained from Channels 1 
(E 0.28 Mev) and 2 (E s 1.2 Mev). This is consistent with the derivation 
of the decay lifetimes obtained by Beall et al .. (1967) and reflects the 
uncertainty in the channel 3 data due to background proton contamination. 

At the higher 1 values (L > 1.8) in this region, the data were Adjusted 
to fit reported measurements I* Van Allen . 1966; Mcllvaln . 1966; Brown . 1966]. 

Finally in region III, which includes all L- space beyond L~ 2.3> . 
the decay lifetime was taken to be only a function of energy. The t'b 
in this latter region arc typically of the order of one month or less and, 
as stated previously, will be neglected in this treatment# Figure 1 Out- 
lines the domains of the regions discussed above for electrons with 
E 2r 0.28 Mev, and indicates the dependence of t on the variables 
appropriate to each domain. 
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In formulating the problem, a four dimensional space is considered in 
which the decay lifetime t is a continuous function of the variables B, L 
and energy E. The boundaries beyond which computation of t is dis- 
continued are: 

(i) The equatorial envelope. 

(ii) The radiation belt terminus. 

(lii) The t ■ 0 surface 

where the equatorial envelope is as defined by Mcllwain (1961); and where the 
radiation belt terminus is designated as the constant electron intensity 
contour of one particle per square centimeter per second as given for trapped 
electrons with energies E * 0.5 Mev by Vette et al .. 1966. 

It is assumed that the terminus is time and energy independent in 3 , L 
space, since any slight changes in its position due to energy dependence or 
secular variations do not affect the derivations and conclusions presented 
in this treatment. The four dimensional space is then subdivided into volume 
elements such that within each volume, the parameter t is a linear function 
of all its variables (i.e., all partia 1 derivatives of t arc constant). 

Since t depends mainly on the magnetic shell parameter L, appropriate 
increments in L have to be taken sr* as to preserve the linearity requirement 
within each volume element. Thus, region I is subdivided into 11 L-intcrvalc 
and region II in 19 L-intervals. Figure 2 shows the dependence of t 
on the variables in region I 


for the three cnerQr channcla of satellite 19&3 38c. In the next tvo 
figures, several cross sections of the t domain in reclon X ore plotted. 
Figure 3 ohowa four 12- L croco cectiono for E 2 0.25, 0.5, 1, and 3 Mev; 
and Figure 4 shove ccveral D-E crone cectiono at discrete L-valucs. In 
the latter figure, the B dependence of t at these low L values is 
clearly evident. Finally, the *r dependence on L and E in region II 
is shown in Figure 5* The slowest decay in this model occurs at L~ 1.6 
for E~ 0.3 Mcv and gradually drifts toward L ~ 1.4 at the higher integral 
energies. Curves for t are not plotted for E < 0.2 Mcv because there is 
evidence that the low energy component (E ~ 0.1 Mev) of the integral 
electron flux does not manifest a perceptible decay r Pflt2cr , 1968; 

Pfltzcr and Wlnckler , 1968]. 


Data Processing 

A generalized treatment for all three variables, such as pertaining 
to region I, will be given here. Hie same treatment, appropriately modified 
to deal with only two variables/ "applies to region II. In general, 

t = t(B;L,E) (1) 

The following relations then hold within each volume element: • 


T« K.(6>,L}E. + AC'&jL) 


K: 

3F- 




M a a C< L + S 



yl + 


(2) 


and 


h ■ r(Ob* er(u) 
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(3) 


Combi nine the expressions ve obtain a general function for the mean lifetime 
vithln the volume element 

T(B,L,e)« «BLE.+pBL-f¥l-E.+ SBEt6.B+^L+^E + 95 (4) 

The coefficients of the variables, B,L and E form a matrix arid arc found 

u 

experimentally, for a particular volume element, byji linear fit to the 
1963 3 OC data. For those intermediate volume elements for vhich coefficients 
cannot be determined due to sparse or insufficient data points, the known 
coefficient values from adjacent bracketing volume elements ore linearly 
extrapolated, so as to match all slopes at the respective boundaries. 

Thus, a continuous functional relationship of the decay lifetime t on its 
variables B,L and E is generated for the inner zone. Table 1 gives the 
values of the 30 sets of coefficients in the respective L~ ranges. 


C»>tvl Ian:*- 


In using the results presented In thio treatment, the following 
cone idem tiono should be taken into account: 

(a) Decay lifetimes should not bo extrapolated backward too 
close to the epoch of the nuclear event because the initial decay (firat 
few weeks) is considerably faster than at .later times, which arc considered 
in the present analysis. 

(b) The estimated maximum value of the error for the calculated 

t is ~ 15 days at intermediate L values (1*3 < L < 2*0). In this region 
the decay lifetime is of the order of several hundred days. At low and 
high L values (L < i .15 and L > 2.2), where the t’s are typically of 

t 

the order of a few weeks, the estimated error is ^ 5 days. 

(c) The functional dependence of t on energy has been derived 
using integral energy channels, and may differ from decay constants 
obtained from differential energy measurements. 

(d) At low energies (K < 0,15 Mev), the'validity of the energy 
dependence of t becomes marginal. There is evidence that electrons of 
these energies do not manifest a decay r Pfltzcr , 1968; Pfltaor and 
Wincklcr , 1968 ] but are mainly part of the natural component of the 
radiation belt. At the higher integral energies, the treatment should 
not be extrapolated beyond E 2 6 Mev* 
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Figure Captions 


Figure 1* Demainn of functional dependence of the decay lifetime 
on B, L and E for E a 0.28 Mev. 

Figure 2. Decay lifetimes in region I (1.15 £ I* $ 1.4) showing the 
B, L dependence of t for the 1963 3&C satellite energy 
channels . 

Figure 3* B-L cross sections of the t domain in region X (1.15 5 
L£ 1.4) for several integral energies. 

Figure 4. B-E crocs sections of the t domain in region X (1.15 £ 
L £ 1.4) at several L values. 

Figure 5* Decoy lifetimes in region XX (1.4 £ L £ 2. 3) for Integral 
energies above 0.2 Mev. 
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